Photon-axion conversion in Active Galactic Nuclei? 
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Axion-Like Particles (ALPs) are the focus of intense current research. We analyze photon- ALP 
conversion in the context of relativistic jet models of Active Galactic Nuclei (AGN) for more than 
100 sources. Contrary to previous claims, we find that this process cannot occur above 100 GeV 
regardless of the actual AGN model and the values of ALP parameters. This result rules out a pro- 
posed strategy to bypass the cosmic opacity above 100 GeV, as apparently required by observations. 
We also show that for some AGN an observable effect can show up in the X and soft 7-ray bands. 
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Introduction - Phenomenological as well as conceptual 
arguments entail that the Standard Model (SM) should 
be regarded as the low-energy manifestation of a more 
fundamental unified theory of all interactions. Invari- 
ably, a theory of this sort involves a set of new particles. 
Therefore, the SM lagrangian is expected to be modified 
by small terms describing interactions among known and 
new particles. Many extensions of the SM which have at- 
tracted considerable interest in the last few years predict 
the existence of Axion-Like Particles (ALPs). They are 
spin-zero light bosons defined by the low-energy effective 
lagrangian 
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where F^" is the electromagnetic field strength, F^i, is 
its dual, a denotes the ALP field whereas m stands for 
the ALP mass. According to this view, it is assumed 
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250 GeV. On the other hand, it is sup- 
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posed that m <^ Gp ~ 250 GeV. Besides than in four- 
dimensional models [J, 0, [3, 0] , ALPs naturally arise in 
the context of compactified Kaluza-Klein theories [5|] as 
well as in superstring theories @, Moreover, it has 
been argued that an ALP with mass m ~ 10"^'^ eV is 
a good candidate for the quintessential dark energy 
The standard Axion [^, [l3| is the archetype of ALPs and 
is characterized by a specific relation between M and m, 
while in the case of generic ALPs M and m are assumed 
unrelated. So, the peculiar feature of ALPs is the tri- 
linear 7-7-a vertex described by the last term in /Ialp, 
whereby one ALP couples to two photons. 

Because of such a vertex, ALPs can be emitted by 
various astronomical objects, and this fact yields strong 
bounds: M > 0.86 • 10^° GeV for m < 0.02 eV ^ and 
M > 10" GeV for m < 10-i°eV Moreover, 



the same 7-7-a vertex makes the interaction eigenstates 
differ from the propagation eigenstates, and so photon- 
ALP oscillations show up when an external magnetic 



(electric) field is present [15|, |l6|, ll7| . 

Several attempts have recently addressed the astro- 
physical implications of ALPs. In this Letter, particular 
attention will be paid to very-high-energy (VHE) effects, 
which can be summarized as follows. Below 100 GeV 
photon propagation throughout cosmic distances is basi- 
cally unaffected by any cosmic background radiation [l8| , 
whereas above 100 GeV Extragalactic Background Light 
(EBL) makes the horizon of the observable Universe 
rapidly shrink as the energy further increases, since the 
process 77 — > e'^e^ becomes an important source of 
opacity [l9|. As stressed above, photon- ALP oscillations 
require the presence of a magnetic field, and in fact spe- 
cific situations have been envisaged, like a magnetic field 
in the source, cosmic magnetic fields [13, [lH along the 
line of sight and the Galactic magnetic field. Actually, 
photon- ALP oscillations have been proposed as a means 
to circumvent the cosmic opacity - as apparently required 
by observations - by relaying upon the fact that ALPs do 
not suffer EBL absorption. Two specific scenarios have 
been put forward: 

• Photon-ALP oscillations are supposed to occur 
only in intergalactic space (22. 23|: 



• Photo-to-axion conversion is assumed to take place 
in the source and to be followed by back-conversion 
in the Galaxy 



Photon-ALP conversion in the source has been envis- 
aged to take place inside the relativistic jets of Active 
Galactic Nuclei (AGN) Yet, no extensive quan- 

titative analysis has been carried out so far. Our aim is 
to fill this gap. We begin by framing our analysis within 
leptonic AGN models, but later we will consider hadronic 
AGN models as well. 

Photon-ALP oscillations - For our needs, it is suffi- 
cient to consider the propagation of a photon beam in 
a constant magnetic field B. In such a situation, the 
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probability that a photon converts into an ALP after a 
distance x can be computed exactly and reads 17 1 



(2) 



with 
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where 
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stands for the plasma frequency of the medium, me is the 
electron mass, is the corresponding electron density, 
Bt is the component of B transverse to the beam and 
— 4.41 ■ 10^'^ G is the critical magnetic field. 

AGN leptonic models - Among the various kinds of 
AGN, blazars are the ones that look brighter in the VHE 
band. They are modelled as a pair of oppositely directed 
relativistic jets - one pointing towards the observer - as- 
sociated with a rotating supermassive black hole. The 
jets are supposed to be powered either by magnetic fields 
twisted by differential rotation of the accretion disk sur- 
rounding the black hole [2§| or by the energy extracted 
from the black hole ergosphere [29|. Photon emission is 
assumed to arise from very energetic electrons contained 
in one or more blobs with typical size of about 10^^ cm 
moving along the jet (30j . In order to explain the strong 
observed variability of jet emission more than one blob 
is needed in some cases. 

A crucial ingredient of this picture is the presence of a 
magnetic field B in the jets, which consists of one or more 
domains larger than the blob [sij. As a consequence, the 
electrons in the blobs radiate through the Synchro Self 
Compton (SSC) mechanism. Namely, relativistic elec- 
trons emit synchrotron radiation up to the X-ray band, 
whose energy is boosted up to the TeV region by In- 
verse Compton scattering on their parent electrons. It 
can also happen that low-energy photons coming from 
outside the blob are upscattered in the same fashion, a 
situation called External Compton (EC) scattering. This 
simple model is able to successfully account for the spec- 
tra of more than 100 blazars. 

We begin by considering the emission from a single 
blob. We use as input parameters in Pj^a{x) the values 
listed in Table 1 of ref. [s^ and we assume that photon- 
ALP conversion occurs inside one single blob (this restric- 
tion will subsequently be relaxed). We first compute the 
survival probability P-y^^ = 1 — P^^a for a photon inside 
a blob and next we convolve it with the pure SSC spec- 
trum, thereby predicting the AGN spectra in the pres- 
ence of photon- ALP oscillations (we are disregarding for 




FIG. 1: We plot By„a for 2251 + 158 vs. E (eV). This figure 
is obtained for M = lO" GeV, m = 10"^° eV. 
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FIG. 2: We plot vF^ (ergcm"^ s" 



3. E (eV) for 2251-^158. 
The solid line is the spectrum as computed in ref. [s^l while 
the dotted line is the spectrum including ALP-Photon oscil- 
lation. Dashed line represents the difference between the two. 
This figure is obtained for M = lO" GeV, m = 10"^° eV. 



simplicity the existence of extragalactic magnetic fields, 
which can give a further contribution and should be con- 
sidered in a more complete analysis pHl. 

We have analyzed more than 100 sources, assuming 
10-1° eV < m < 10-6 eV and M = 10" GeV as illus- 



trative values. We find P-y^a > 0.1 for 12 of them and 
P^^a > 0.5 only for 2 sources: 0528-1-134 and 2251-^158. 
Let us focus for definiteness on 2251 + 158. We get a non- 
negligible conversion only within a finite energy range 
and the corresponding behaviour of P-y^a is plotted in 
Fig. IT]). Surprisingly, the conversion happens to take 
place at energies much lower than the ones previously re- 
ported in the literature 24, 2^ 2^, 2^. Moreover, while 
the existence of the low-energy cut-off is commonly ap- 
preciated, the relevance of the high-energy one is often 
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neglected in practice [S^l- Still, the existence of such a 
high-energy cut-ofF can be very important in some as- 
trophyisical settings. Indeed, in the present discussion 
the above behaviour of P^^a emerges for all considered 
sources, and in no case does a nonvanishing P-y^a show 
up for E > lOGeV. 

Remarkably enough, the shape of P^^a exhibited in 
Fig. ([T]) can be understood in an intuitive fashion. To 
see this, we note that from Eq. ^ it follows that a 
necessary condition for P^^^a to be maximal is 
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In order to work out the relevance of Eq. ([5]) , it proves 
instrumental to define 
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so that we have m, ~ m for m ^ Wpi and to, ~ Wpi for 
TO <C ujpi. In addition, we set 
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Correspondingly, Eq. ^ can be rewritten as 
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where e{m — tOpi) denotes the sign of its argument. Now, 
by combining Eqs. (O and ([9]), it is straightforward 
to check that P-y^a is maximal in the so-called strong- 
mixing regime defined hy E^, < E < i?,,, which is indeed 
the range A — C shown in Fig. We stress that E^, 

and i?,, are independent of the size of the magnetized 
region in which photon-ALP takes place, and moreover 
-B*, is also independent of to and uj-p\: these facts play a 
key-role throughout the subsequent investigations. 

We proceed by extending our analysis to the many- 
blob case. What happens can be sketched as follows. 
Focusing one one generic blob, on top of the single-blob 
behaviour discussed above, the EC emission from blobs 
closer to the AGN centre has to be taken into account. 
Furthermore, we use the same assuniptions and scaling 
laws for the blobs as outlined in ref. [33|. We have han- 
dled this more complicated situation by means of a nu- 
merical code [3J|. In this case, P-y^a turns out to be 
vanishingly small even for E^ < E < i?,,, simply be- 
cause the decrease of the blob size - needed to keep the 
total resulting luminosity constant - drastically reduces 
photon-ALP conversion according to Eq. 
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FIG. 3: We plot Pj^a for E = lOMeV photons in a B (G) 
vs. linear dimension (cm) plane. White corresponds to no 
conversion at all while black corresponds to P-y^a > 0.9. This 
figure is obtained for M = lO" GeV, m = 10"* eV, Ue = 
10^ cm"^. 



AGN hadronic models - A natural question is whether 
Pj^a can be sizeable for E > 100 GeV within AGN 
hadronic models. Remarkably, this issue can be settled 
by merely performing an analysis of the ALP parameters 
space. The result is exhibited in Fig. We see that 
P-ywa is non-vanishing inside a series of elongated sets 
which move in the parameter space towards wider and 
less magnetized AGN jet regions as E increases. Specifi- 
cally, in order to have a sizeable P-y^a for E > 100 GeV, 
Eq. ^ requires 



Bt < 4.5 • 10" 



10^0 GeV\ 
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Taking as before M ~ 10" GeV, we find that Bt < 
4.5 • 10~^ G is needed. Yet, hadronic models demand 
lOG < S < lO^G inside the jet [H, thereby implying 
that Py^a is always negligibly small for E > 100 GeV. 

Discussion - So far, we have supposed that photon- 
ALP conversion takes place inside the blobs, but this 
might not be the case. Hence, we have run the multi- 
blub simulation for 10 blobs allowing for the conversion to 
happen also in the intra-blob region of the jet. The result 
is depicted in Fig. (|4]). As it is evident, now Pj^a can be 
sizeable and quantitatively the situation happens to be 
(just by chance) practically the same as in the single-blob 
case. It goes without saying that the same argument can 
be applied to the single-blob case, possibly resulting in 
an enhanced P.y^a [35!]. 

On the basis of our previous finding, one might think 
that a magnetic field smaller than Bt ^ 4.5 • lO^^G 
could give rise to a successful conversion if the model- 
dependent constraint WG < B < 10^ G were relaxed. 
However, this turns out to be an impossible task. For, it 
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both in the source and in intergalactic space, as stated in 
the text (such an analysis is currently being performed 
for the X and soft 7-ray emission from the sources listed 
in Table 1 of ref. 
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FIG. 4: We plot uF^ (erg cm"^ s"^) vs. E (eV) for a ran- 
domly generated source with 10 emitting blobs. The solid 
line represents the pure SSC spectrum, while the dotted line 
gives the same spectrum including photon- ALP oscillations. 



does not matter to have > 100 GeV if P^^a is small 
by its own. As it is clear from Eq. we also have to 
require the sine argument to be great enough so as not 
to kill the conversion probability. It is easy to see that in 
the strong mixing regime the sine argument is Bx/2M, 
which is of order 1 provided that B is coherent on a 
scale of about 10^^ cm. Yet, observations [31] tell us that 
at distances greater than lO^'* cm from the black hole, 
B becomes chaotic and looses its coherence, thereby de- 
stroying photon- ALP oscillations. For the same reason, 
we conclude that no conversion for E > 100 Gev takes 
place in the radio lobes opening at the end of the jet. 

Conclusions - We have analyzed photon- ALP conver- 
sion in the context of AGN relativistic jet models. We 



have found that - contrary to previous claims |24l . l25l . 
- this process cannot occur above 100 GeV regardless of 
the actual AGN model and the values of ALP param- 
eters. Thus, only photon-ALP oscillations occurring in 
intergalactic space can successfully circumvent the cosmic 
opacity in the VHE band 13, 23]. We have also shown 
that for some sources an observable effect can show up 
in the X and soft 7-ray bands [s^. 

Note added - Just before the submission of this Letter, 
we have become aware of a very recent paper [37| which 
discusses photon-ALP conversion in the source and in in- 
tergalactic space. This paper considers only 2 sources: 
3C279 and PKS 2155-304 (see theh Table 2). In either 
case, they find that for E > 100 GeV no conversion inside 
the source takes place, in agreement with our findings. 
So, the problem of the cosmic opacity for these sources 
can indeed be side-stepped only within the scenario ad- 
vocated in refs. [H, |2^. For E < 100 GeV things are 
different, and we agree with ref. 37[ that a full-fledged 



analysis should take into account photon-ALP conversion 
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